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A measurement of direct photons in p+p collisions at ^/s — 200 GeV is presented. A photon 
excess above background from ty° — > 7 + 7, r\ — *• 7 + 7 and other decays is observed in the transverse 
momentum range 5.5 < pt < 7 GeV/c. The result is compared to a next-to-leading-order pertur- 
bative QCD calculation. Within errors, good agreement is found between the QCD calculation and 
the measured result. 

PACS numbers: 13.85.Qk,25.75.Dw 



Measurements of particle production at large trans- 
verse momenta (pt) in hadronic interactions provide the 



* Deceased 

tPHENIX Spokesperson zajc@nevis.columbia.edu 



possibility to test perturbative quantum chromodynam- 
ics (pQCD). Neutral-pion production in p+p collisions 
at y/s = 200 GeV in the range 2 < p T < 13 GeV/c 
measured by the PHENIX experiment at RHIC can be 
well described by next-to-leading-order (NLO) pQCD pj . 
This comparison, however, relies on the choice of the 
parton-to-pion fragmentation function. Measurement of 
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direct-photon production provides a more direct test of 
pQCD. Quark-antiquark annihilation (qq — > jg) and 
quark-gluon Compton scattering (qg — > (77) contribute 
to direct-photon production at leading order pj. Due 
to the latter process, which dominates the production, 
the measurement of direct photons can be used to ob- 
tain information on the parton distribution function of 
the gluon inside the proton. 

Previous direct-photon measurements in p+p collisions 
were made up to energies of s/s = 63 GeV (see e.g. 
jiJ [3 0) For p + p collisions direct-photon data are 
available at considerably higher en ergies , y/s — 546, 630 
GeV HUH up to y's = 1800 GeV jalllj. At these ener- 
gies NLO pQCD calculations describe the direct-photon 
data within about 20%, although systematic differences 
in the spectral shapes were observed 0]. At energies 
below y/s = 63 GeV the agreement between NLO pQCD 
and data is generally worse. With phenomenological ap- 
proaches based on soft-gluon radiation of the incoming 
parton, which leads to an additional transverse momen- 
tum fcx, a better description of the data can be obtained 
The Relativistic Heavy Ion Collider (RHIC) at the 
Brookhaven National Laboratory provides p+p collisions 
at energies between the existing data sets, allowing better 
constraints on the processes affecting incoming partons. 
A further incentive to study direct-photon production in 
p+p at yfs — 200 GeV comes from the measurement of 
direct photons in collisions of gold nuclei at the same 
center-of-mass energy per nucleon-nucleon pair ll l| . 
In central Au+Au collisions high-px neutral-pion pro- 
duction is suppressed 0] which is due to energy loss of 
the scattered quarks and gluons in the hot and dense 
fireball created in these collisions (jet quenching) [l6| . 
As direct photons are not subject to the strong interac- 
tion they should not be suppressed in the jet-quenching 
model. In this context, the p+p direct-photon results 
serve as a baseline against which possible nuclear effects 
can be identified. 

The data presented in this report were collected during 
the 2001-2002 run period (Run 2) of RHIC. The neutral- 
pion cross sections obtained from this data set were pub- 
lished in . The unpolarized neutral-pion spectrum and 
the unpolarized direct-photon spectrum presented here 
were obtained by averaging over proton bunches with 
varying vertical polarization delivered by RHIC. 

Direct photons and background photons from decays 
7T° — > 7 + 7 and r\ — » 7 + 7 were measured with the 
electroma gne tic calorimeter (EMCal) of the PHENIX ex- 
periment |17|. The background from charged particles 
was subtracted with the aid of a layer of multi-wire pro- 
portional chambers with pad readout (PC3) which was 
located directly in front of the EMCal. The minimum- 
bias trigger was provided by two beam-beam counters 
(BBC) which were also used to determine the collision 
vertex. In addition to the minimum-bias trigger condi- 
tions, a high-px photon trigger was used, derived from 
the analog energy signal measured with the EMCal. 

The two BBC's were located symmetrically around the 



nominal interaction point at ± 1.44 m along the beam- 
line. The BBC's subtended the pseudorapidity range 
± (3.1 — 3.9) with full azimuthal coverage. The colli- 
sion vertex was determined by measuring the difference 
of particle arrival times in the two BBC's. This analy- 
sis was restricted to events with a vertex in the range 
± 30 cm. The BBC's were calibrated as a luminosity de- 
tector with absolute luminosity measurements based on 
the van der Meer scan technique 18]. With these scans 
the cross section for firing the BBC minimum-bias trigger 
was determined to be 21.8+ 2.1 mb p|. Thus, roughly 
50% of the inelastic p+p events satisfy the minimum- 
bias trigger condition if an inelastic p+p cross section of 
42 mb at y/s = 200 GeV is assumed. 

The PHENIX EMCal comprises two arms each with 
4 sectors [l^. The EMCal consists of two different sub- 
detectors, a lead-scintillator calorimeter (PbSc, 6 sectors) 
and a lead-glass calorimeter (PbGl, 2 sectors). Each sec- 
tor covers a pseudorapidity range of \rj\ < 0.35 and an 
azimuthal range of A0 w 22.5°. Each PbSc (PbGl) sec- 
tor is highly segmented and consists of 72 x 36 (96 x 48) 
individual detector modules, called towers, with a lateral 
size of 5.5 x 5.5 cm 2 (4x4 cm 2 ). With a radial distance 
of the sectors to the beamline of roughly 5 m this corre- 
sponds to a segmentation of A<j> x A77 as 0.01 x 0.01 such 
that the two decay photons of a ir° are well separated 
up to neutral-pion momenta of px ~ 20 GeV/c. The 
different detection mechanisms of the two sub-detectors 
(measurement of scintillation light in PbSc and detection 
of Cherenkov photons in PbGl) result in a different re- 
sponse to hadrons. Thus, the PbSc and PbGl provide 
photon measurements with different systematic uncer- 
tainties. The energy calibration of the detector was ob- 
tained from the position of the 7r° invariant-mass peaks. 
A ~ 4% (~ 5%) shift of the 7r° peak position due to 
energy smearing in conjunction with the influence of the 
steeply falling ir° px spectrum was taken into account in 
the PbSc (PbGl) calibration. The calibration was cor- 
roborated by correlating the EMCal energy with the mo- 
mentum of electrons measured with the PHENIX track- 
ing detectors and, in case of the PbSc, by measuring the 
energy deposited by minimum-ionizing particles. From 
these studies the systematic uncertainty of the energy 
measurement was estimated to be less than 1.5%. In a 
direct-photon analysis it is essential to exclude bad de- 
tector modules ("hot towers") which might give rise to 
spurious direct-photon signals. A detailed quality assess- 
ment was carried out to identify such towers. 

The EMCal high-px trigger (called 2x2) was based 
on the analog energy signal measured in 2 x 2 groups of 
adjacent EMCal towers (called trigger tiles). The av- 
erage threshold of the trigger corresponded to an en- 
ergy signal of 0.75 GeV. The probability as a function 
of the photon px to fire the trigger was determined by 
Monte Carlo simulations which included the variation of 
the trigger tile thresholds, the EMCal detector response, 
and the geometry of the active trigger tiles. This trigger 
efficiency was confirmed with minimum-bias data. The 
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photon trigger efficiencies for PbSc and PbGl reached a 
plateau above pt = 1.5 — 2 GeV/c at the limit of about 
0.78 expected from the number of active towers and 2x2 
trigger tiles. The high-px-trigger photon sample was used 
above pt = 3 GeV/c in the final spectrum. 

Another EMCal trigger which did not require a coin- 
cidence with the minimum-bias trigger was used to ac- 
count for the bias on the particle measurement due to the 
minimum-bias event selection. To this end the fraction of 
ir°'s measured with this EMCal trigger for events which 
in addition satisfied the minimum-bias trigger condition 
was determined to be / = 0.75 ± 0.02. The unbiased 
photon and neutral pions cross sections were then deter- 
mined by dividing the total number of measured photons 
and neutral pions by this number. 

The minimum-bias data sample in this analysis con- 
sisted of 16.7 million events, corresponding to an inte- 
grated luminosity of 0.77 nb _1 . About 1 in 47 minimum- 
bias events also satisfied the 2x2 high-p T trigger condi- 
tion. The 18.7 million analyzed 2x2 events thus corre- 
sponded to an integrated luminosity of 40.3 rib -1 . 

The first step in the direct-photon analysis was to de- 
fine a sample of direct-photon-candidate hits. An EM- 
Cal hit was rejected as a direct-photon candidate if it 
formed an invariant mass in the ir° or r\ range with other 
hits in the same or adjacent sectors. The invariant-mass 
window was 110 < m 77 < 170 MeV/c 2 for the n° and 
500 < to 77 < 620 MeV/c 2 for the r], corresponding 
roughly to a ±2ct window around the observed ir° and r\ 
peaks. To keep the rate of accidental rejections of genuine 
direct photons low it was required that the partner hits 
had a transverse momentum of pt > 0.4 GeV/c. This cut 
effectively corresponded to a px-dependent upper limit 
on the energy asymmetry a = \E 7 i — E^/iJS^i + E^) 
in the rejection procedure. In spite of this requirement, 
a small fraction of genuine direct photons is rejected. 
This was studied by inserting artificially generated direct- 
photon hits into real events. In order to keep the hit 
multiplicity constant, a randomly selected real hit was 
removed from an event in this procedure. It was found 
that the loss of genuine direct photons was less than 2% 
for pt > 3 GeV/c. The final direct-photon spectrum 
was corrected for this effect. In order to increase the 
chances of finding the partner photon for a ir° or r\ de- 
cay photon, direct-photon candidates were required to lie 
within a restricted fiducial area which was defined by a 
minimum distance of 16 (20) towers to the edge of the 
detector for the PbSc (PbGl). For example, for 7r°'s with 
Pt = 4 GeV/c and the requirement that one decay pho- 
ton has a pt > 0.4 GeV/c the average distance of the 
decay photons in tower units is ~ 8 (~ 11) for PbSc 
(PbGl). With the chosen fiducial area basically all de- 
cay photons from neutral pions with pt ^ 4 GeV/c can 
be tagged. Monte Carlo studies showed that the rejec- 
tion of direct-photon candidates based on the tt° and r\ 
tagging lead to a reduction of background photons from 
hadron decays in the fiducial area of about a factor of 
2 for pt > 5 GeV/c. Some direct-photon analyses only 



measure isolated direct photons for which the total trans- 
verse energy or the number of charged tracks in a cone 
centered around the direct photon is required to lie below 
a threshold. No such cut was used in this analysis. 

In order to reduce the background from hadronic hits 
in the EMCal, cuts were applied on the lateral shower 
shape and on the time-of-flight of the hits. The remain- 
ing contamination of charged particles was subtracted 
on a statistical basis by employing the PC3 as a charged- 
particle veto detector. The intrinsic efficiency of the PC3 
for detecting a charged particle was higher than 99% 
and the active PC3 area in the EMCal acceptance was 
roughly 90%. PC3 hits were projected onto the EMCal 
surface using a straight line given by the PC3 hit and the 
event vertex. An EMCal hit within a certain veto radius 
was counted as a charged hit. The chosen veto radius 
decreased with increasing pt and for pt > 0.8 GeV/c a 
constant value of 15 cm was used. The fraction of charged 
hits was corrected for random associations with the help 
of a mixed-event technique. The charged-particle back- 
ground in the direct-photon-candidate sample was ~ 15% 
around pt = 5 GeV/c for both PbSc and PbGl. A large 
fraction of these background hits, however, comes from 
photon conversion in the field-free region between the 
vertex and PC3. The photon loss due to conversion was 
calculated based on the material budget up to PC3. The 
photon conversion probability was 4.1% for the 2 PbSc 
sectors in the East Arm of the central spectrometer, 5.3% 
for the 4 sectors in the West Arm, and 7.4% for the PbGl. 
These conversion losses were taken into account in the 
final photon cross section. The correction for the con- 
tamination of the raw spectrum of neutral EMCal hits 
with neutrons and anti-neutrons was determined with a 
detailed GEANT simulation [2jJ. In the case of the PbGl 
calorimeter the simulation was based on the creation of 
Cherenkov photons in order to achieve a realistic descrip- 
tion of the detector response. The background from neu- 
tral particles was found to decrease with pt and was al- 
ready less than 1% for pt > 2 GeV/c for both PbSc and 
PbGl. 

The geometric acceptance and the efficiency of the pho- 
ton detection were calculated with a Monte Carlo simu- 
lation. The efficiency takes the distortion of the direct- 
photon-candidate spectrum due to energy smearing into 
account. Moreover, it corrects the small (< 5%) loss of 
photons due to the shower shape and time-of-flight cuts. 
The final direct-photon spectrum was corrected for the 
difference between the average direct-photon cross sec- 
tion within a finite pt bin and the value of the cross 
section at the bin center. 

With the described corrections the unbiased differ- 
ential cross section J^f^. t = Ed 3 a/d 3 p for the direct- 
photon candidates calculated from the minimum bias 
data sample reads 



j3 -i -i / ' s~i s~i jucana 

0.(7 1 1 t-^rcco ' L'conv ' ^loss direct 7 

d 3 p £ 2-kpt f Apx • Ay 



5 













■ (data) 








; + 


Ybckg ' (calculated background) 




r * 








r * 






: , , , i , 


+ 


- — _i 






) 2 


4 6 8 


10 


12 



p T (GeV/c) 



FIG. 1: (color online) Ratio of direct-photon candidates to 
the 7T° spectrum. The histogram represents the expected 
background signal from the Monte Carlo calculation. The 
error bars represent the statistical error and the boxes the 
systematic error. 



TABLE I: Systematic uncertainties of the neutral-pion spec- 
trum, the direct-photon-candidate spectrum, and the mea- 
sured and simulated "f/n° ratios at px = 6.75 GeV/c. 



7T° error source 


PbGl 




PbSc 


Yield extraction 


5% 




5% 


Yield correction 


8% 




6% 


Energy scale 


9% 




8% 


Total 


13% 




12% 


Tdfrcct error source 


PbGl 




PbSc 


Non-7 background correction 


4% 




4% 


Yield correction 


6% 




5% 


Energy scale 


9% 




9% 


Total 


12% 




11% 


Total error 


PbGl+PbSc combined 


7direct/ 7r 




9% 




7bckg/7T° 




4% 





where ^direct 7 i s the total number of direct-photon 
candidates in a px bin Apx and rapidity bin Ay; C rcco 
is the acceptance and efficiency correction for photons; 
Cconv is the correction for photon conversions; Ci OS s is 
the correction for the loss of genuine direct photons in 
the 7T° and r\ tagging; / = 0.75±0.02 is the fraction of the 
unbiased direct photon yield which is measured under the 
minimum bias trigger condition; and C is the integrated 
luminosity for the analyzed data sample. The high-px 
triggered sample required an additional correction for the 
efficiency of this trigger for photon detection. 

The pt spectrum of the direct-photon candidates con- 
tains direct photons as well as remaining background 
photons from hadron decays. These background pho- 
tons mostly come from ir° and 77 decays for which one 



decay photon misses the detector. At a representative 
bin of pt = 6.75 GeV/c about 93% of the background 
photons originate from ir° and r\ decays, the remaining 
background photon come from decays of other hadrons 
like lo and rj' . The background was calculated with the 
same Monte Carlo code that was used for the acceptance 
and efficiency calculation. The Monte Carlo code took a 
parameterization of the measured tt° spectrum as input. 
The spectra of r\ mesons and other hadrons with photon 
decay branches were assumed to have the same sha pe as 
the 7T° spectrum as a function of wit = \/v\ + m o ( m T 
scaling 0,|2^]). The rj/n invariant cross section ratio 
as a function of itit was taken as 0.48 ±0.1 which was 
confirmed by the measured 77 spectrum. 

The dominant systematic uncertainty of the neutral- 
pion spectrum and the direct-photon-candidate spectrum 
came from the uncertainty of the energy scale and the un- 
certainty of the yield correction. At px = 6.75 GeV/c, 
the 1.5% uncertainty of the energy scale resulted in a 
9% uncertainty in the photon yield. The yield correction 
included the correction for energy and position smear- 
ing of the detector, for photon losses due to particle- 
identification cuts, for photon conversions, and for the 
detector acceptance. For the neutral pions an additional 
5% uncertainty came from the extraction of the 7r°-peak 
content. In the case of the photon measurement the 
uncertainty due to the subtraction of charged and neu- 
tral backgrounds was taken into account. In the ratio 
7dircct/ 7r ° °f tne direct-photon-candidate spectrum and 
the neutral-pion spectrum systematic uncertainties par- 
tially cancel. Monte Carlo studies showed that the uncer- 
tainty of this ratio at px = 6.75 GeV/c due to a possible 
non-linearity of the energy scale was ~2%. 

For the determination of the direct-photon spectrum 
the expected background photons from hadronic decays 
need to be subtracted from the spectrum of direct- 
photon candidates. To this end the ratio i? 7 = 
(7dfroct/ 7r °)/(7bckg/7'' ) of the measured direct-photon 
candidates to the calculated background was determined. 
The direct-photon spectrum was then calculated as 

7direct = (1 - R- 1 ) ■ 7^ (2) 

The relative systematic uncertainty of the direct- 
photon cross section was calculated as the quadratic 
sum of the relative uncertainties of the two factors in 
Eq. [21 The factor 1 — i?" 1 contains the (statistical 
and systematic) significance of the direct photon signal. 
When multiplying with J^f^. t , on iy the systematic un- 
certainties that cancelled in the ratio i? 7 are added {e.g. 
the energy scale error). The overall normalization un- 
certainty from the luminosity determination was 9.6%. 
The estimated systematic uncertainties for the measured 
neutral-pion spectrum and the measured direct-photon- 
candidate spectrum are shown in Table[l]for a representa- 
tive bin (pr = 6.75 GeV/c) of the spectrum. The system- 
atic uncertainty of the direct-photon measurement was 
corroborated by comparing results obtained for the dif- 
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FIG. 2: (color online) Measured cross section and NLO 
pQCD calculations for direct-photon production in p+p colli- 
sions at y/s —200 GeV. The normalization error of 9.6% is not 
shown. The two data points plotted with an arrow indicate 
the beginning of the low- and high-pT ranges where the direct 
photon signal is consistent with zero. The upper edges of the 
arrows indicate an upper limit (90% confidence level) for the 
direct photon cross section calculated from the statistical and 
systematic uncertainty. 



TABLE II: Invariant differential cross section for direct- 
photon production in p+p collisions at y/s = 200 GeV. Asym- 
metric uncertainties (cri ow , (Thigh) are given for the cross sec- 
tion. The absolute normalization error of 9.6% is not in- 
cluded. 



PT 


£d 3 (r/d 3 p 


stat. 


error 


sys. 


error 


(GeV/c) 


(mb GeV _2 c 3 ) 


O"low 


Ohigh 


O"low 


Thigh 


5.75 


5.6L10" 7 


50% 


42% 


53% 


54% 


6.25 


2.68-10" 7 


75% 


56% 


55% 


56% 


6.75 


2.37- 10" 7 


59% 


42% 


37% 


38% 



ferent photon-identification criteria. Moreover, the indi- 
vidual PbGl and PbSc results were found to agree within 
systematic errors. 

To make the best use of the available Run-2 statis- 
tics the photon and w° spectra from PbSc and PbGl 
were combined for the final result. The ratio of the 
acceptance- and efficiency-corrected direct-photon can- 
didate spectrum to the measured ir° spectrum is de- 
picted in Fig. ^ In addition, the Monte Carlo calculation 
for the ratio of the expected background photons to the 
7T° spectrum is shown. At high px the statistical signif- 
icance of the direct-photon-candidate spectrum is weak. 
However, around pt ~ 6— 7 GeV/c there is clear evidence 
of a photon signal above the background. 

The extracted invariant direct-photon spectrum is 



shown in Fig. [5] and the numerical values are given in 
Table im The experimental result was compared to NLO 
pQCD calculations El El El El El. 

which used the 

CTEQ6 parton distribution functions |29| and the GRV 
parton-to-photon fragmentation function |3C| . There are 
in general two mechanisms for the production of direct 
photons: the direct contribution from elementary scat- 
tering processes of quarks and gluons, described in the 
introduction, and the contribution from photons which 
are produced in the fragmentation of quark or gluon jets. 
The latter is a long-distance process which is not pertur- 
batively calculable. It is described by a parton-to-photon 
fragmentation function which is determined experimen- 
tally. Since no isolation cut was used in the data anal- 
ysis the pQCD calculation in Fig. [21 includes contribu- 
tions from the direct production mechanism and the frag- 
mentation mechanism. The separation of short-distance 
and long-distance processes in the pQCD calculation in- 
troduces unphysical renormalization, factorization, and 
fragmentation scales. Identical values for all three scales 
were used in the pQCD calculation. In Fig. [3 results are 
shown for three choices of the scales (/i = pt, [i = pt/2, 
and n = 2pt). The theoretical and experimental results 
agree within the large uncertainties of the data points. 

Prior to the Run-2 p+p beamtime PHENIX took data 
from Au+Au collisions at ^/snn = 200 GeV. A clear 
direct-photon signal was observed in mid-central and cen- 
tral Au+Au reactions The strong suppression of 
7r°'s and rj's in Au+Au significantly reduced the number 
of background photons and eased the extraction of the 
direct-photon signal. The p+p NLO pQCD was used as 
a baseline reference for the interpretation of the Au+Au 
result. In contrast to neutral pions no sign of a suppres- 
sion of direct photons in Au+Au collisions was found. 
The direct-photon measurement presented in this paper 
supports the use of the NLO pQCD calculation as a ref- 
erence for the results measured in Au+Au. 

In summary, a small but significant direct-photon sig- 
nal has been observed at mid-rapidity in p+p collisions 
at y/s = 200 GeV. The measured direct-photon cross sec- 
tion is in agreement with pQCD calculation, albeit within 
large errors. 
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